Dynamics of a tunneling magnetic impurity: Kondo effect induced incoherence 
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We study how the formation of the Kondo compensation cloud influences the dynamical properties 
of a magnetic impurity that tunnels between two positions in a metal. The Kondo effect dynamically 
generates a strong tunneling impurity-conduction electron coupling, changes the temperature de- 
pendence of the tunneling rate, and may ultimately result in the destruction of the coherent motion 
of the particle at zero temperature. We find an interesting two-channel Kondo fixed point as well 
for a vanishing overlap between the electronic states that screen the magnetic impurity. We propose 
a number of systems where the predicted features could be observed. 
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Tunneling of a heavy particle or some collective degree 
of freedom in a dissipative environment has been the sub- 
ject of intense theoretical and experimental research in 
the past and is by now reasonably well-understood . 
Maybe the most intriguing case is that of ohmic dissipa- 
tion, where the particle usually couples to electron-hole 
excitations of a metallic environment. In this case the 
bare tunneling amplitude Ao of the particle is strongly 
renormalized due to the dissipative environment, and be- 
comes temperature dependent. In the simplest scenario, 
where the tunneling occurs between two sites, the ef- 
fective tunneling displays a power-law behavior over a 
wide range of temperatures, A(T) ~ A (T/u ) a , with 
ujq a high-energy cut-off of the order of the Debye fre- 
quency, and a a dimensionless coupling constant KB. 
The renormalization of the tunneling amplitude is a con- 
sequence of Anderson's orthogonality catastrophe || : At 
any position the presence of the particle generates a 
screening cloud that consists of an infinite number of 
electron-hole excitations. The formation of this huge 
'electronic polaron-cloud' slows down the particle, in- 
creases its mass, and thus decreases its tunneling am- 
plitude. Depending on the specific value of the coupling 
a, the dynamics of the particle can be of three differ- 
ent kind: (a) If the coupling is small the particle moves 
with slightly damped coherent oscillations between the 
two positions at T = 0. (b) For 1/2 < a < 1 the motion 
of the particle becomes incoherent, while for even larger 
couplings (c) the particle becomes localized and cannot 
move from one well to the other. 

At low temperatures and larger values of a, more com- 
plicated dynamical processes involving simultaneous tun- 
neling and electron scattering (electron assisted tunnel- 
ing) may become more relevant and can, in principle, 
delocalize the particle or lead to other types of strongly 
correlated states @JnJ. In this paper, however, we focus 
on the rather generic case of 'slow' TLS's ||, where the 
energy scale generated by these processes is extremely 
small, and they can therefore be neglected. 

In the present paper we focus our attention to the 



very interesting but poorly understood case of a tunnel- 
ing magnetic impurity coupled to an ohmic environment 
jll],[L2|]. Possible examples of such a system include a 
magnetic impurity tunneling between an STM tip and a 
metallic surface jl^| , a Kondo impurity in an amorphous 
region a spin 1/2 quantum soliton interacting with a 
metallic environment pJl, or charge tunneling in double 
quantum-dot systems [|16| , [L7| . In all these cases the spin 
of a magnetic impurity couples to the conduction elec- 
trons through an exchange interaction. This exchange 
interaction becomes renormalized as the temperature is 
lowered and may lead to the dynamical formation of a 
Kondo compensation cloud. Since the tunneling particle 
has to drag this cloud with itself, the exchange coupling 
will lead to a dynamically generated temperature depen- 
dent renormalization of the tunneling amplitude. In the 
present paper we study this interplay between the Kondo 
effect and the orbital motion of the TLS in detail. 

For the sake of simplicity we focus our attention to 
the simplest possible case of a TLS, where tunneling 
takes place between two positions only, R±. Further- 
more, though we also discuss the role of asymmetry to 
some extent, we mostly focus on spatially symmetrical 
TLS's. We show, in particular, that the Kondo effect as- 
sociated with the magnetic degrees of freedom leads to 
a strong temperature- dependence of the exponent a, and 
may eventually induce a decoherent state. Destroying 
the Kondo cloud with a magnetic field one could then be 
able to drive the particle back to the coherent regime. 
The orbital motion, on the other hand, may lead to the 
appearance of a two-channel Kondo state under special 
conditions, where the impurity tunnels very fast back 
and forth and forms a Kondo state with the conduction 
electrons at both positions. 

We describe the TLS by the tunneling Hamiltonian 



ff tun = -(A T a! +A,T*) 



(1) 



where the two pseudospin states T z = ±1/2 correspond 
to the two tunneling positions, Ao is the tunneling matrix 
element, and A z describes the asymmetry of the TLS. 
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As a first approximation, we consider the exchange in- 
teraction between the tunneling particle and the electrons 
as completely local: 

H int = ]T J g P g S ■ ( 2 ) 

Here P± = 1/2(1 ± 2T Z ) projects out the TLS states 
at positions R±, J± is the exchange coupling at these 
positions, S denotes the spin operator of the impurity, 
and a stands for the Pauli matrices. The field opera- 
tors, = / e 4kR ' ± c ] J c/j d 3 k/(27r) 3 , create conduction 
electrons at R± with spin \i. Including more complicated 
non-local interactions in Eq. (||) such as spin-assisted tun- 
neling has not changed our results [ |lT|]l6f| . 

First we consider the case of a symmetric TLS with 
J_ = J j. — J and A 2 = 0. Following the standard pro- 
cedure fllCyTil , the relevant conduction electron degrees 
of freedom can be represented simply by one-dimensional 
fermion fields , c pall , obeying canonical anticommuta- 
tion relations, {cp afi , c],,^,} = 2ir 5(p-pO<W<W, and 
described by 

H ei = [ ^ VFP 4<v W ■ ( 3 ) 

In Eq. (||) the radial momentum p is measured from the 
Fermi momentum hp, vp is the Fermi velocity, a — {e, o} 
is the parity, and u denotes the spin. Introducing the 
fields = J_% p Cpa^dp/V^, Eq. (§) becomes 

H int = |s(l + F)*tff* e + |£(1 - F)*t^ 

+ gSVl - F 2 T z (tf „ + . (4) 

Here F = sm(kpd)/kpd measures the overlap of the 
states V&i, with d = |R+ — R_| the tunneling distance, 
and g — Jk F /2ir 2 . For d = this Hamiltonian obvi- 
ously reduces to the single-channel Kondo model. In the 
following we set Vp = kp = K = 1. 

To study the zero temperature dynamical properties 
of the tunneling Kondo spin we used Wilson's numerical 
renormalization group (NRG) Jh|. In this very accurate 
numerical technique one constructs a series of Hamilto- 
nians, -ff/v, which are diagonalized iteratively. Having 
obtained the many-body eigenstates and energies of Hn 
one can use them to calculate physical quantities at an 
energy scale T,w ~ uj n ~ A~ (JV+1)/2 , with A w 3 a 
discretization parameter. Our results were obtained by 
keeping the lowest 250 states in each iteration. To obtain 
accurate results we exploited the following symmetries: 
(i) parity (ii) global spin rotations, and (iii) a hidden 
SU(2) symmetry, related to electron-hole symmetry [ pp[ . 

We computed the T = impurity spin and pscudospin 
spectral functions, g l = — Imlx^w)}, (with XoC^) the 
Fourier transform of the retarded response function) , us- 
ing their Lehmann representations: 



{?o(w > 0) = Y Kn|O*|0)| a 5(B n -Ea- to). 

n 

Here | 0) (| n)) is the ground (nth excited) state of the 
system with energy Eq (E n ), and O l denotes the opera- 
tors T 1 and S\ 
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FIG. 1. Logarithmic plot of various spectral functions dis- 
cussed in the text for A = 3, g = 0.144, and A = 2.31 ■ 1CT 5 . 
The energy scales Tk and A are also indicated. 

Our results are summarized in Fig. |l|. Let us first fo- 
cus on the somewhat peculiar case of F = 0, where \E' + 
and do not overlap. In this case we can observe two 
distinct cross-overs: The first takes place at the Kondo 
energy Tk ~ e^ 1 / 29 and corresponds to the formation of 
a local Kondo state at the actual position of the TLS. 
Above Tk all spectral functions behave as 1/u), indicat- 
ing that all correlation functions are constant for times 
shorter than 1/Tk- Below Tk the spin spectral func- 
tion becomes linear. Performing a Hilbert transform on 
can see that this corresponds to a constant impurity sus- 
ceptibility ~ 1/Tk- The z-component of the pseudospin 
spectral function is practically unaffected by the forma- 
tion of the Kondo compensation cloud: at this time scale 
tunneling events are very rare, and the particle can be 
considered as immobile. This interpretation is confirmed 
by a detailed analysis of the finite size spectrum |n| . The 
logarithmic slope of the x-component of the pseudospin 
spectral function, however, does become modified: this 
change is related to the dynamical renormalization of the 
tunneling amplitude by the formation of the Kondo com- 
pensation cloud. 

Though A <C Tk in our computations, A is a rele- 
vant perturbation, and leads to a second cross-over at a 
renormalized value of the tunneling amplitude, A* , where 
the TLS freezes into the even tunneling state. For F = 
a two-channel Kondo state is formed below A*, as con- 
firmed by the analysis of the finite size spectrum [116). 
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This is most easily understood by observing that the last 
term of Eq. (g) flips the TLS between the T x = ±1/2 
states, and can therefore be dropped below A*. Here 
.Hint becomes simply the two-channel Kondo Hamilto- 
nian, and a two-channel Kondo state is formed in the 
spin sector The spin spectral function and the z- 
component of the pseudospin spectral function become 
constant below A*, implying the logarithmic divergence 
of the spin- and T 2 -susceptibility at low temperatures, 
XS (T) ~ A*ln(T K /T)/T%, and xf,(T) ~ ln(T K /T) / A* , 
and the asymptotic behavior of the time-dependent cor- 
relation functions (S l (t)SH0)} ~ (T z (t)T z {0)) ~ 1/t. 

These results immediately imply that the external 
magnetic field and asymmetry are both relevant opera- 
tors of scaling dimension 1/2 at this two-channel Kondo 
fixed point |21jj . The two-channel Kondo state is thus 
extremely unstable. It crosses over to a Fermi liquid 
state also for any finite overlap F: For small values of 
F a third crossover occurs at an energy T* well below 
A*, however, for generic F's this second crossover takes 
place almost simultaneously with the crossover at A* and 
only a small kink remains from the two-channel Kondo 
behavior of F = 0. 

To determine the scale A* and the temperature depen- 
dence of the tunneling rate we invoked renormalization 
group arguments. The Hamiltonian involves three di- 
mensionless parameters: g, F, and Aq/ujo, with ojq the 
high-energy cutoff. Well above A*, A /ujo is small and 
the effective tunneling at energy scale uj or temperature 
T ~ uj, A(w), satisfies the following scaling equation: 



dlnA(uj) 
d\n(uJo / uj) 



-a{g(w),F), (wo»A(w)). (5) 



Above A* the impurity is immobile. Therefore the effec- 
tive Kondo coupling g(uj) does not depend on the other 
two parameters, and is a universal function of uj/Tk- 
In this regime F can be completely transformed out of 
the Hamiltonian and is therefore constant. Thus a — 
ci(uj/Tk,F) is a universal function of lo/Tk and F for 
u)q 3> A. This function is related to the time-dependent 
correlation function as (T x (i)T x (0)) - l/t 2a( -^ (t ~ 
1/lo) ^lj. Assuming that the retarded and the time or- 
dered correlation functions have the same singular be- 
havior, this immediately implies that, within logarith- 
mic accuracy, a can be determined from the logarithmic 
derivative of Qrp as 



1 ,dhx q% 

a(uj) Ri -(— 

y ' 2 V d\nuj 



+ 1), (w»A*). 



(6) 



We plotted this universal function in Fig. For large 
frequencies a ~ 0, meaning that the tunneling ampli- 
tude remains unrenormalized above Tk- Below Tpc, on 
the other hand, it scales to an overlap-dependent con- 
stant, a < (F). As shown in the inset, a< coincides with 
the Anderson orthogonality exponent K for a maximally 
strong scatterer with a phase shift = 6i = n/2 22 1, 



1 tan S a VT^F^ \ 2 

atan — , ] . (7) 



V 7 ! + F 2 tan 2 <5 CT 



Thus, far below Tjc (but still above A*) the Kondo im- 
purity acts as a maximally strong potential scatterer in 
agreement with Nozieres' Fermi liquid picture pgfl. 
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FIG. 2. Top: Energy- and overlap-dependence of the 
anomalous dimension a of the dimensionless tunneling am- 
plitude for F = 0, 0.1, 0.2 0.3, 0.4, 0.5, and 0.6 (top to 
bottom), as determined from the logarithmic derivative of 
the spectral function g%. Bottom: Temperature dependence 
of the normalized tunneling rate, R(T) = R(T)/Ag. F de- 
creases from top to bottom. We used A = 3, g = 0.144, and 
A = 2.31 • 10" 5 in both figures. 

A quantity of primary interest is the renormalized tem- 
perature dependent tunneling rate, which can be exper- 
imentally determined by performing real-time measure- 
ments 0. By dimensional analysis, for an ohmic heat 
bath it is given by R(T) - A(u = T) 2 /T §. We com- 
puted it by integrating Eq. (|^) (Fig. bottom). As the 
most striking consequence of the Kondo effect, the loga- 
rithmic slope of R{T) changes at T w Tk- 

The energy scale A* is determined by the condition 
A(w = A*) w A*, leading to 



A (C 



A \ Q </( 1_Q <) 
Tk~) 



(8) 



with C a constant of the order of unity. The constant 
a< also characterizes the dissipative non-equilibrium dy- 
namics of the TLS below Tk 0- In Fig. | we show 
the rescaled spectral function, A* 2 q^(uj) / 'uj , related to 
the real part of the retarded response function. With- 
out dissipation, a < — (F = 1), the tunneling of the 
TLS is entirely coherent: the TLS oscillates between the 
two positions without damping, and the spectral function 
consists of two Dirac delta's. For F < 1 the coherence 
peak broadens: The oscillations become exponentially 
damped and at very long time scales (at T = 0) the corre- 
lation function behaves as (T z (t)T z (0)) ~ 1/t 2 . For even 
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smaller values of overlap the peak becomes completely 
invisible, and the formation of the Kondo compensation 
cloud suppresses the coherent motion. 
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FIG. 3. Rescaled spectral functions and the evolution of 
the coherence peak with increasing overlap. We computed A* 
from Eq. (0) with C = 0.0916 and T K = e" 
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In the general case, the TLS model is more complex 
than the one we discussed until now. The TLS is not 
symmetrical: A 2 ^ 0, and the couplings J± in Eq. (||) 
are not equal. The difference between J± leads to two 
subsequent Kondo effects at T^'s associated with the two 
positions, and consequently changes in the logarithmic 
slope of R(T) twice, while a finite A 2 generates a new 
energy scale, below which the TLS freezes into one of the 
states T z = ±1/2 ||). 

In Eq. (^|) we only took into account the exchange in- 
teraction. In reality, the TLS and the electrons interact 
through a local potential scattering as well. However, 
for a Kondo impurity this potential scattering is rela- 
tively small (of the order of the exchange coupling) in 
the exchange scattering-channel and since it remains es- 
sentially unrenormalized, its effect can be neglected com- 
pared to that of the exchange interaction. Potential scat- 
tering in other scattering channels may, however, be still 
present and shift a by a temperature-independent value 
a(T) — > a(T)+a > . Therefore, for small overlaps one can 
easily be in a situation, where the Kondo effect drives the 
TLS from a coherent ground state to an incoherent state 
with a > 1/2 (see Fig. |). 

For small enough Kondo temperatures, the value of 5 a 
in Eq. um can be changed by an external magnetic field 
p3| , which gradually destroys the Kondo effect, thereby 
changes the logarithmic slope of R(T) to its value above 
Tk, and restores the coherent TLS motion p6| . 

In the experimental realizations mentioned in the in- 
troduction it is the Kondo spin that moves and carries 
the compensation cloud. Our discussion can be, however, 
easily generalized to the case of a TLS that happens to be 
close to a Kondo impurity |24|,^5| . In this case Friedel os- 
cillations generated by the TLS modify the local density 
of states at the Kondo impurity and thus the exchange 
coupling will depend on the position of the TLS, leading 
to a Hamiltonian similar to Eq. (^|). 

It is also an interesting question whether the two- 
channel Kondo behavior found can possibly be observed 



in realistic systems. Clearly, one needs an almost per- 
fectly symmetrical TLS and a very small overlap F. In 
most cases F is typically large, excepting the magnetic 
impurity tunneling between an STM tip and a surface, 
where the electronic states coupled to the spin hardly 
overlap at the two positions and F « 0. In order to ob- 
serve the two-channel Kondo behavior, one has to mini- 
mize the TLS asymmetry too: The asymmetry A z can be 
tuned to zero with the external voltage, while J + w J_ 
can be assured by using sufficiently rounded tip that is 
made of the same metal as the surface. 
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